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Abstract. The new representation of the Coulomb three-body wave function via 
the well-known solutions of the separable Coulomb two-centre problem 4>j{£,T)) = 
Xj(£)Yj(r]) is obtained, where X,(£) and Yj(rj) are the Coulomb spheroidal functions. 
Its distinguishing characteristic is the coordination with the boundary conditions of the 
scattering problem below the three-particle breakup. That is, the wave function of the 
scattering particles in any open channel is the asymptotics of the single, corresponding 
to that channel, term of the expansion suggested. The effect is achieved due to the 
new relation between three internal coordinates of the three-body system and the 
parameters of <frj(£,r]). It ensures the orthogonality of 4>j(£,n) on the sphere of a 
constant hyperradius, p = const, in place of the surface R = \x% — x\\ = const 
appearing in the traditional Born-Oppenheimer approach. The independent variables £ 
and 77 are the orthogonal coordinates on that sphere with three poles in the coalescence 
points. They are connected with the elliptic coordinates on the plane by means of the 
stereographic projection. For the total angular momentum J > the products of 
<j>j and the Wigner D-functions form the hyperspherical Coulomb spheroidal (HSCS) 
basis on the five-dimensional hypersphere, p being a parameter. The system of the 
differential equations and the boundary conditions for the radial functions ff(p), the 
coefficients of the HSCS decomposition of the three-body wave function, are presented. 
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1. Introduction 

The complete set of the solutions of the Coulomb two-centre (CTC) problem is widely 
used as a basis for the representation of the Coulomb three-body (CTB) wave function 
for a long time (pQ and references therein). One of the most attractive features of that 
basis is its simplicity. First, due to the high symmetry of the Coulomb field, the CTC 
problem admits the separation of variables, and the basis elements are the completely 
factorized functions. Second, the co-factors composing the basis elements, so-called 
Coulomb spheroidal functions (CSF), are the well investigated functions related to the 
Heun class. The principal results on CSF and the references are presented in [2]-[l]- 

The traditional way to use the CTC basis for the analysis of the CTB system with 
positive and negative charges {Z X Z 2 > 0, ZiZ 3 < 0) is the adiabatic representation [1] 
which is known also as the Born-Oppenheimer (BO) approach, the perturbed stationary 
states method, the molecular orbitals method. In that approach the distance R between 
the particles 1 and 2 is considered as an adiabatic parameter, and the basis functions 
are orthogonal at fixed R. The numerous helpful properties of CSF allow to construct 
the effective numerical algorithms for the calculation of bound states and cross sections 
of elastic and inelastic collisions for various CTB systems. The method is successfully 
used not only for the case of heavy particles 1,2 and a light particle 3 p], [5]- [9], but as 
well for light particles 1,2 (electrons) and a heavy particle 3 (nucleus) [T0]-[T3"]. 

Nevertheless, the BO representation in the traditional form (i.e. with R as 
a parameter) has a well-known imperfection which decreases the efficiency of the 
calculations of the scattering processes. The point is that the BO expansion is not 
coordinated with the physical boundary conditions of the scattering problem. It means 
that at R — > oo the wave function of the scattering particles in any open channel is 
formed as the sum of the infinite number of components, the selected terms of the BO 
series calculated at large R. That sum includes, in particular, the basis functions tending 
at large R to the atomic wave functions of the closed channels. Thus, the terms of the 
BO expansion, generally speaking, are not divided into the groups corresponding to 
open and closed channels. The peculiarity pointed out becomes apparent in the system 
of radial equations: the matrix elements coupling different basis states may remain 
nonzero at large R. That leads to the complication of the boundary conditions for the 
radial functions at R — > oo [8] and to the slow convergence of the expansion. The 
incoordination of the adiabatic representation with the physical boundary conditions 
of the three-body (not necessary Coulomb) scattering problem is peculiar to both the 
quantum approach and the semiclassical one. To overcome the difficulties in the frames 
of the semiclassical treatment the inclusion of "electron translation factors" is used 
[H], [15]. The various methods of the improvement of the adiabatic representation for 
the scattering processes are developed also in the quantum approach ([Kj, p2] and 
references therein). 

In contrast to the traditional BO representation the adiabatic hyperspherical (AHS) 
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one [IB]- [22] is coordinated with the boundary conditions of the three-body scattering 
problem. In that approach the hyperradius p is used as an adiabatic parameter. The 
wave function of any open channel is the asymptotics of the single, corresponding to 
that channel, term of the expansion. The AHS approach is successfully applied to the 
calculations of the bound states, the cross sections of the elastic and inelastic collisions 
and the resonances in various CTB systems [18], [20]-[28j. 

However, the AHS basis elements are the essentially more complicated functions 
than the BO ones. In the general case they are the non-factorized functions of five 
independent variables. In addition, for the Coulomb three-body systems the AHS energy 
terms as functions of p have the numerous avoiding crossings which make difficulties for 
the numerical calculations [20], [22]. These avoiding crossings are closely connected with 
the exact crossings of the CTC adiabatic terms [22] . Thus, in the AHS approach the 
high symmetry of the Coulomb field, instead of the simplification, paradoxically leads 
to the additional difficulties. The avoiding crossings are removed in the original version 
of the hyperspherical approach developed in papers [30]- [34]. In these papers the basis 
completely factorizing in the special four-pole elliptic coordinates on the hypersphere is 
suggested. The co-factors composing the basis elements in that approach are essentially 
more complicated functions as compared with CSF, they relate to the Heun class only 
at zero energy [33] . 

The goal of the present paper is to obtain the representation of the CTB wave 
function which, on the one hand, contains only the well-known CSF and, on the 
other hand, is coordinated with the boundary conditions of the scattering problem at 
energies below the threshold of the three-particle breakup. In other words, to obtain the 
representation combining the advantages of the AHS approach (the coordination with 
the boundary conditions) and the advantages connected with the specific character of 
the Coulomb field (the complete factorization of the basis and the simplicity of its 
components, CSF). That task is interesting, firstly, for the development of the efficient 
computing methods and, secondly, for the completeness of the theoretical description 
of CSF: it is desirable to present various ways to use these well-known functions in the 
CTB problem. 

The above goal is achieved in the paper with the help of the new relation between 
the parameters of CSF and the internal coordinates of the CTB system. The suggested 
CTC basis functions, just as the AHS basis functions, form the orthogonal set on the 
sphere p = const in the three-dimensional space of the internal variables, in place of the 
surface R = const appearing in the traditional BO approach. 

The difference of our expansion from the traditional BO one can be outlined briefly 
for the simplest case of zero total angular momentum, J = 0, in the following way. Let 
Xi, x 2 , and x 3 be the position vectors of the particles in the centre-of-mass frame 

miXi + m 2 x 2 + m 3 x 3 = 0. (1) 

In the case J = the three-body wave function \I/ depends on three independent internal 
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variables, for example r l5 r 2 , R: 

ri = \x 3 - xx\, r 2 = \x 3 - x 2 \, R=\x 2 -x 1 \. (2) 

In the traditional BO approach three independent variables are, in fact, R and the 
prolate spheroidal coordinates £ G [1, oo) and 77 G [—1,1] defined by the relations 

£ = (n + r 2 )/R, v = (r 1 -r 2 )/R. (3) 

The variable R is considered as a parameter of the basis, and the decomposition of the 
three-body wave function has the form 

* = J2fAR)UR\^v), (4) 

j 

where the basis functions (j>j(R\£,r)) are the factorized solutions of the CTC problem 
with the intercentre distance R. These solutions as the functions of £ and 77 form the 
complete orthogonal system at fixed R. 

In our approach the independent variables are £, 77 and the hyperradius 

p = \jHm\x\ + 777,2X2 + ^3) 

= \J%m^m^r\ + m 2 m 3 r\ + m\m 2 R?)j(m\ + m 2 + 777,3), pG[0, 00). (5) 

It is considered as the parameter of the basis. The decomposition of the three-body 
wave function has the form 

* = c(p,z, v )Y,fAp)Mp\^v), (6) 

j 

where C(p, £, 77) is the special weight factor. Basis functions 4>j{p\£, 77) are the factorized 
solutions of the CTC problem with the modified charges and the intercentre distance 
proportional to the hyperradius p. They form the complete orthogonal system at fixed 
p and coincide (up to a constant factor) with (f>j(R\£,r)) at large R. 

In the case J > the basis suggested consists of the productions of CSF and 
the symmetrized Wigner functions. We call it "the hyperspherical Coulomb spheroidal 
(HSCS) basis". The formulae are more complicated but the principal idea is the same. 

The principal notations and the starting equations are presented in section 2. 
Section 3 is devoted to the construction of the HSCS basis. In section 4 the definition 
of ^-matrix in the appropriate representation is presented, and the asymptotical 
expressions for the radial functions in terms of the matrix elements of S'-matrix are 
obtained. These formulae are deduced from the properties of the HSCS basis and the 
general relations defining ^-matrix, without any analysis of the system of the radial 
equations. That system and the statement of the scattering problem at energies below 
the three-particle breakup are presented in section 5. 
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2. Starting equations 

We consider the system of three particles with charges Z\ > 0, Z 2 > 0, Z 3 = — 1 
and masses mi, m 2 , m 3 . The case of the identical particles 1 and 2 is not considered. 
The atomic units are used. The Jacobi coordinates R, r are expressed in terms of the 
position vectors of the particles x%, x 2 , a? 3 (CQ) by the formulae 

m x x 1 + m 2 x 2 

R = x 2 -x 1 , r = x z ■ , (7) 

mi + m 2 

where M and p are the reduced masses: 

m 1 m 2 m 3 (mi + m 2 ) 

iW ; > f 1 ; ; • l°J 

mi + m 2 m i + m 2 + TTI3 

The Jacobi coordinates r a ,R a , which are suitable at large R for the clusterization a 
(a = 1 corresponds to the atom (1,3) plus the distant particle 2; a = 2 corresponds to 
(2,3)+l), are defined by 

m a x a + m 3 x 3 



x 3 — Rat — \ 1 



,a-l 



(9) 



ra a m 3 „, m 3 _ Q m a + m 3 

A*a= 1 , M a = ; ; , a = 1,2. (10) 

m a + m 3 mi + m 2 + m 3 

The factor (— l) a_1 (jUJ) ensures the identical directions of R a and at large R. 
The Hamiltonian of the system in the coordinates R, r has the form 

H = T + V, (11) 

T = —A R - —A r , (12) 

2M H 2/i ' v ; 

= ^2 _ Z i Z * f 13 ) 

i2 |r + RM/mi| |r-HM/m 2 |' v ; 

In the hyperspherical approach six variables 

$G[0,27r), 0e[O,7r], y?G[0,2vr), pG[0,oo), xe[0,7r], E [0, vr] (14) 

are generally used. Here {$, B, 0} are the Euler angles of the rotating frame with the 
third axis directed along R. The variables $ and <p are the spherical angles of r in the 
rotating frame. The hyperradius p and the hyperangle x are expressed in terms of R 
and r as 

p 2 = 2(MR 2 + pr 2 ), pe[0,oo); (15) 

tan X /2 = (p/M^r/R, X ^ M- (16) 

The angles <£>, O, <p are the external coordinates defining the orientation of the triangle 
formed by three particles. The variables p, \i $ are the internal coordinates of the 
system. In the hyperspherical coordinates the kinetic energy (Tl2|) takes the form 



1 9 5 <9 4 fd , d , 2 \ J 2 -2(1, J) 

p° ap ap p z sm xvax c^X / p cos z x/^ 
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The operators J and I in (j!7p are the total angular momentum and the angular 
momentum of the third particle respectively. The expressions for I 2 , J 2 , and J ■ I 
are given by 



l 2 = 
J 2 = 
J I 



1 d gin ^_g 
sin i? d$ d-d 



1 



2 



sin 1? <9<£> 2 



+ 



sin <p- 



sin 6 96 S 09 

sin^ 9 
— ^-(cos6 
sm 6 

<9 2 



COS d 2 



sin 2 6 V<9$ 2 

9 . d 



d 2 d 2 n „ 
+ ~ — ~ - 2 cos 6 



d 2 



+ cos <f cot 6 



dp 2 
d_ 

dip 2 



<H>d(f) ' 



(18) 
(19) 



+ 



cot i? 



<9<£> 2 



(20) 



dQdp sin 6 d$d(p 

In what follows we use instead of x the variable 

t = tanx/2, ie[0,oo). (21) 

6ne can consider t and •& as the polar coordinates on the half-plane presenting the 
stereographic projection of the hemisphere with the radius 1, the latitude \ and the 
longitude $ on the equatorial plane (figure 1). Three independent variables t, i? and p 
are the spherical coordinates of the vector t in the rotating frame. That vector differs 
from the relative vector r/R [35] by the mass factor: 

t = tr/r = (n/M) 1/2 r/R. (22) 

For the formulation of our approach it is suitable to write the three-body 
Hamiltonian H (flTl) . ([121) . (Ti~3l in terms of the variables $, 9, tp, p, t, Taking into 
account the expression for the six dimensional elementary volume 

dRdr = (AM/j,)- 3/2 gt 2 sin ■& sin 6dpdtdtfd^d6d$, (23) 

S = p 5 (l + t 2 )- 3 , (24) 

one can present the kinetic energy T (j!2p as 

T = q-^Tq^ 2 , (25) 



-1/2/- 

<9p 2 



+ 



1 + t 2 )A t (l + t 2 ) + (1 + t 2 ){J 2 - 2J ■ I) + 



where 



A, 



19 J 
-t 2 



t 2 df dt t 2 ' 
Using (|22|) and the notations 

t a = (— l) a t a R/ R, 



(/xM) 



1/2 m -1 



a 



1,2, 



(26) 



(27) 



(28) 
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we obtain V ( 1131) in the form 

:i+t 2 ) i/2 



V 



\t-t r 



p 



(2M) l/2 Z 1 Z 2 - 



t 2 + t 2 - 2(-l) a tt a cos $ 



t — til 

1/2 



|t-t 2 | 

a = 1,2. 



(29) 



(30) 



The three-body wave function ^ = ^(p, £, i?, $, O, <p) satisfies the Schrodinger 
equation 

= EV. (31) 

Using the transformation (1251) one can present it in the form 

M = EV, (32) 

where the transformed wave function \1/ and the transformed three-body Hamiltonian 
H are given by the relations 

* = g l /2* = p B ' 2 (l + 1*)- 3 ' 2 % (33) 



H = g 1/2 Hg- 1 ' 2 



d 2 



dp* 



'i + t 2 



(j 2 - 2 J ■ Z) + 



4p 2 



(34) 



-p- 2 (l + t 2 )A t (l + t 2 ) + V. (35) 

The transformed AHS Hamiltonian h (|35|) acts on functions of t, <p and depends 
on p parametrically. Its eigenfunctions are used as a basis for the representation of the 
three-body wave function in the reduced version of the AHS approach [22], [23]. (In the 
general version [19], [20] the AHS Hamiltonian /i ad includes two last terms of (26) and 
acts on functions of five variables.) 



3. The hyperspherical Coulomb spheroidal (HSCS) basis 

3.1. The appropriate Coulomb two-centre (CTC) Hamiltonian in t-space 

It is seen that /i ad ([35]) can be considered as the Hamiltonian of the particle with the 
variable mass m(t) moving in the field V in t-space: 

h ad = -(2m)- 1 / 2 A t (2m)- 1 / 2 + V, 2m(t) = p 2 (l + t 2 )~ 2 . (36) 

Here the form of the kinetic energy ensures the self-conjugacy. The pure discrete 
character of the spectrum of /i ad is explained by the fast decrease of the mass m(t) 
at large t, rh(t) ~ t -4 . Every eigenfunction of h ad coincides at large p with the wave 
function of certain bound state of the atom (1,3) or (2,3) [22]. 

For the generating of the basis functions in our approach we also use the 
Hamiltonian of the particle moving in t-space, but the mass of the particle is fixed. 
Namely, we use the traditional CTC Hamiltonian 
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It describes the motion of the particle with the mass 1/2 and the charge -1 in the field of 
two effective charges Z\ and Z 2 located in the points t\ and t 2 of t-space at the distance 

|t 2 =h+t 2 = (p/M) 1/2 (38) 

from each other. In our approach the effective charges Zi = Zi(p) and Z 2 = Z 2 (p) are 
chosen so that every eigenfunction of the discrete spectrum of h (13"?]) at large p coincides, 
as well as in the case of /i ad , with the wave function of some bound state of the atom 
(1,3) or (2,3). The simplest choice is the linear function Zj(p) for which we easy obtain 

Z a = pZ a {2p) l ' 2 {l + t 2 a y 3/2 = pZ a 2 l ' 2 p^ 2 p-\ a = 1, 2. (39) 

Indeed, under that condition each eigenfunction of the discrete spectrum of h ( 1371) 
at p — > oo is localized (with the exception of the case Z\ = Z 2 , see below) in the 
neighborhood of one of the centres and coincides with the eigenfunction of the one- 
centre Hamiltonian [3]. For example, in the asymptotic region 

R^oo, \x 3 -x 1 \=r 1 = 0(l), (40) 
or, in the hyperspherical variables, 

p -> oo, \t - ti| = (p/MY^n/R = n(2//xi) 1 /V/p + 0(/T 2 ), (41) 
the corresponding one-centre Hamiltonian is given by 

- 2) ■ (42) 

ii n J 

The expression in the parenthesis coincides with the Hamiltonian of the atom (1,3). 
Hence, the eigenfunction of h (l37j) localized at p —>■ oo in the domain (|40p . (j4"Ij) coincides 
up to a constant factor with the wave function of some bound state of the atom (1,3). 
The indices of the state are unimportant here. 

The formula for Z 2 is proved in exactly the same way. 

As it is seen from equation (I39p in the case 

Z x /Z 2 = (pz/ptf 2 (43) 

the effective charges are equal, Z\ = Z 2 . In that case the Hamiltonian (1371) 
has an additional symmetry, and its eigenf unctions are divided into symmetric and 
antisymmetric ones. Thus, every eigenfunction is localized at p — ► oo near both centres. 
For the construction of the basis in that case it is suitable to use the corresponding 
linear combinations of symmetric and antisymmetric eigenfunctions which are localized 
at large p in the vicinity of one of the centres. That transformation does not present 
any difficulties. Taking into account the abstract character of the condition (j4"3"j) for the 
non-identical particles 1 and 2 we do not consider that case and suppose hereinafter 

Zi ? Z 2 . (44) 




Hyperspherical Coulomb spheroidal representation in the three-body problem 



9 



3.2. Coulomb spheroidal functions on the sphere p = const 

The equation denning the spectrum of the Hamiltonian h(p\i) ( |37l) . 

(h-e)ip = 0, (45) 

is investigated in details [21 [3]. The problem admits the separation of variables in the 
prolate spheroidal coordinates £ G [1, oo) and rj G [—1, 1], 



V 



-'(I* 



til + I* 



*2|) 



iT 1 (ri + r 2 ) 



iT 1 (n 



(46) 
(47) 



The eigenfunction ipf m (p\t), which corresponds to the eigenvalue e = Ej m < in the 
discrete spectrum and to e > in the continuous spectrum, is presented as the product 

#m(p|t) = <Mp|£, V)e ±imip /V^, (48) 



(49) 



Here m = 0, 1,2, ... is the absolute value of the projection of the total momentum on 
the rotating axis. The multi-index j is defined for the discrete spectrum (e < 0) and for 
the continuous one (e > 0) by the equation 



J 



(50) 



(e,n v ), e>0, 

where and n v are the number of nodes of radial X(£) and angle Y{rf) functions 
respectively. These functions, the so-called Coulomb spheroidal functions (CSF), satisfy 
the system of equations 



A-p 2 (£ 2 -l) + a£ 



— (1 - 7? 2 )— + A - p 2 (l - rf 2 ) + bq 



x(t) = o, 



Y(rj) = 



(51) 
(52) 



di] drj 

and corresponding boundary conditions [3j. Here A is the separation constant. The 
parameters a, b and p are connected with the parameters of the Hamiltonian h (j3"7|) by 
the relations 

P lr7 .3/2 7 3/2, 



^(i 1 + z 2 )(t 1 + t 2 ) 



^2 



VP? 
p 

^2pM 



{Z2P2 2 



Z2P2 

z^T) 



1 



t 9 



(53) 
(54) 
(55) 



— = 1 

2 W ' ^ 2V M' 

CSF defined by ( l5Tl) - lj55l) differ from traditional CSF [3], [1] by the new connection 
of the parameters a and b with the coordinates and the masses of three particles (1531) . 
( }54"|) . The most important point is that our a and b are proportional to p instead of R 
in the traditional approach. 
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The functions <pj m (1491) at fixed m and all possible j form the complete set in the 
space of functions of two internal coordinates {t, $} or {£,??}) at fixed third 

internal coordinate, p = const. The normalization accepted in the paper is defined by 
the equation 

U ^ jKn'M^ £<0 ' , KR , 

< famWj'm >= < (56) 

{(>n v n' v d{£ -£), £>0, 
where < f\g > denotes the two-dimensional integral 

1 Jsm 2 x/2 



cos 4 y/2 
o ' o 



OO 7T / , ^ OO 1 



= J t 2 dt J sin W f*g = filial / de / d7] ^ - rf) f* g . (57) 
00 1 
The completeness condition takes the form 

j 

00/00 °° 

= £ I £ 0n £ n„m(p|^,^)0n £ n„m(p|^,V) + / e n„m(p|£, ^)0en,m(p|^, V)^ 

= 8(t a + t 2 )- 3 (^ 2 - r/T 1 ^ " 0*fa - »/)■ (58) 

It is convenient to consider the variables p, \ and d as the spherical coordinates (the 
radius, the latitude and the longitude respectively) in the three-dimensional space with 
Cartesian coordinates x = p sin x cos??, y = p sin x sin z = pcosx- The net £ = const, 
r] = const on the sphere p = const in that space is presented on figure 1 (that figure 
is taken from the monograph [36J, only the axes with the notations corresponding to 
our problem are added). It can be obtained by the stereographic projection of the 
sphere on the plane z = as a preimage of the traditional orthogonal elliptic net on 
the plane. Three poles correspond to the coalescence points of the three-body system. 
As the stereographic projection is the conformal mapping the net on the sphere is the 
orthogonal one. 

3.3. HCS basis on the five- dimensional hypersphere 

The three-body Hamiltonian (fill) has three commuting integrals of motion pQ: the 
square of the total angular momentum J 2 , its projection on the third axis of the 
laboratory frame J3 and the inversion of all Jacobi coordinates P. It is suitable 
therefore to use for the construction of the basis the symmetrized Wigner D- 
functions D^ n ($, 0, <p) which are the eigenf unctions of these operators. The functions 
D^ x m (^, 6, if) are, in addition, the eigenf unctions of the operator J' 2 where J3 is the 
projection of the total angular momentum on the third axis of the rotating frame: 

J 2 D J K x m = J(J + l)D J K x m , J 3 D J K x m = KD^ m , (59) 
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Figure 1. The three-pole orthogonal net £ = const, r\ = const on the sphere 
p = po = const in the three-dimensional space x = p sin % cost?, y = psinxsini?, 
z = pcosx; y > in our problem; £ and r\ are given by (|46p . (|47|) . The poles 
correspond to the coalescence points of the three-body system. The stereographic 
projection with the centre x = z = ~Po onto the plane z — gives the usual elliptic 
net on that plane. For p = 1 the polar coordinates on the plane z = coincide with 
the independent variables t and $ of the CTC Hamiltonian (|37|) . 



-rl2p.J\ _ m ^nJA pn JA — \ n JA \ ui 

J 3 ^ifm — '' 4 -^i^m' r u Km ~ AU Kmi A ~ 3=1 - 

./A 



,2 n JA 



!JA 



(60) 



The normalized symmetrized D-functions D^ n ($, 0, <£>) is given by the expression [22 



D 



Km 



+ A(-1 



\J+m-Q J 



K-m 



(61) 
(62) 



where 

A J Km = (~l) K [(2 J + 1)/(1 + 5 0m )] 1/2 /(4tt). 
The standard Wigner functions D J Km (Q, 6, <£>) are defined via functions d Km (Q) 

D J Km (^ 0, = exp HK$)4 m (6) exp (-irn^), (63) 

2"7T 7T 27T 

|d$| sinOdG J d v D J K ^D J K ^, = 6jj,5 KK ,6 mm ,6 xy . (64) 



o o 



We define the hyperspherical Coulomb spheroidal (HSCS) basis on the 5- 
dimensional hyper sphere 



n = {t,>&,*,e,<p} = {z,7i,$,e,<p}, (65) 
(t + 1 ) 3 

dn = t 2 sin^dtd^sin6d$ded^ = (^ 2 - r7 2 )d^dr/sin9d<l»d9dy (66) 
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as the system of the common eigenfunctions of the operators J 2 , J3, J'3, P and the 
Coulomb-two-centre Hamiltonian h (1371) depending on the parameter p : 

^\ P \n) = D J K x m ($, e, V ). (67) 

Due to the completeness and the orthonormality of the system {$j„ A } (EZD the 
three-body wave function \1> (1331) with fixed J, K, A is presented in the form of the 
HSCS-decomposition 

y JKX ( P ,n) = P 5 / 2 (i + t 2 )- 3/2 y JKX (p,n) = J2f JKXjm (p)®j* x (p\ty 1 (68) 

jm 

where the sum by j includes both the discrete spectrum and the continuum (see (1581) ). 
The coefficients of the decomposition (168|) . the radial functions /j^f A (p), are given by 
the integral 

f JKX jm{p) = J ^ {p ^JKX {pmdn (69) 

4. S-matrix and radial functions 

4-1. Basis functions of the discrete spectrum at large p 

Every eigenfunction of the discrete spectrum of h ( 1371) at large p tends to the wave 
function of the bound state of the atom (1,3) (a = 1) or the atom (2,3) (a = 2). 
The detailed analysis [3] shows that the limiting atomic states are the states with the 
definite parabolic quantum numbers in the rotating frame. Thus, the basis functions 
( l49l) can be enumerated by the indices a (the number of the nucleus forming the atom), 
n (the principal quantum number), s (the first parabolic quantum number) and m (the 
absolute value of the projection of the atomic orbital momentum on the axis R): 

0n e n,m(p|f , V) = 4>ansm{p\ii V)- (70) 

The formulae determining the one-to-one correspondence 

{rig, n v , m} <-> {a, n, s, m} (71) 

for different cases are presented in the monograph [3]. 

To obtain the factor connecting the asymptotics of 4> a nsm.(p\£,, v) an d the normalized 
atomic function in the region 

p -> oo, r a = 0(1) (72) 
one have to use the asymptotical formula 

I* — *x| = r a (2/p a ) l / 2 p/p + 0(p- 2 ) (73) 
and relations ( |56l) . ( 1571) defining the normalization of basis functions. The result is 

<t>ansm{p% V) = (PaP 2 /^f'^smK, K)- (74) 

p— »oo 
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Here 4>£$ sm {r a , $a) is the normalized atomic function of the atom a without the factor 
exp (±im(p) / v^27r 

'r 2 a dr a J sin^d^l^^L^,^)! 2 = 1, (75) 

o 

d a is the angle between r a and R . The atomic function is factorized in the parabolic 
coordinates, but we conserve the notation 4>anlm( r a^a) for convenience. 

Using equation (17^1) we obtain the asymptotics of the HSCS basis function ^nsm 
(j67p in the form 

KnL(p\ti V,*,Q,<P) = (/i Q p 2 /2/i 2 ) 3/4 ^($,e,^,r a ,^), (76) 

p— >-00 

where 

F Q J ^($, 9, ^ r a , a ) = D J K x m ($, 6, <p)^ sm (r a , d a ) 

= Km [d J - Km (9) exp(-zm^) + A(-l) J+m ^ K _ m (e) exp(zm^)" 

xap^^SKr,,^. (77) 



It is seen that F^ sm is the combination of two atomic functions with opposite 
projections of the orbital momentum on the rotating axis R with coefficients depending 
on $ and 9. Functions F^n Sm corresponding to the same atom (a = a') satisfy the 
orthogonality condition 

2n n 2n oo 2n 

jd$ j sin6d0 J dip J r 2 a dr a J sm M^ a F^F^Z 



= 0~JJ'$KK'&\\'&mm'$nn'$ss>- (78) 

4-2. S-matrix in the (ansm) -representation 

The product of -F^nsm d^Zj) and the ingoing (outgoing) wave 

^™($, ©, V, r a , # a ) • i^exp [±*(A; an i? a - 7cm )] (79) 

presents the state of the colliding atom and the remaining particle ((l,3)+2 for a = 1 
or (2,3)+l for a = 2) at large distance i? a between them. Here k an is the relative 
momentum 



k an = pM a (E - E an ) = M a v an , (80) 
E an is the energy of the isolated atom a, the logarithmic phase 7 Qn is given by 

7cm = 7an(R a ) = (Z a ~ l)Z 3 - a M a k~* \og2k an R a . (81) 

The asymptotics of any three-body wave function at energy below the threshold of the 
three-particle breakup is presented as the combination of the products (1791 . 

Consider the three-body system with fixed J,K,X. These indices are omitted 
hereinafter. According to the general theory [37] for the definition of iS-matrix we 
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have to consider the three-body wave function which contains the single ingoing wave in 
the input channel {a, n, s, m} and the outgoing waves in all open channels {a', n', s', m'} 

exp [-i(k an R a - j an )} 



_ -*- ansm 7-1 j — 



ST da'n's'm'v exp [j(k Q , n ,R a , - 7a , n ,)] 

2^ b ansm * a'n's'm 1 5 T== • \P^) 

tta' V Va'n' 



a'n's'm' 

' ansm 



Here the coefficients S^ s ^ l m are the matrix elements of the operator 

§=SI, (83) 

where S is the S'-matrix, and J is the inversion operator which changes the sign of the 
relative coordinate R a at fixed r a . It acts on any function of 0, (p, d a as follows 

IY ($, 9, <p, # a ) = Y ($ + 7T, IX - 0, -p, 7T $a) • (84) 

The matrix elements of S'-matrix in any representation can be obtained from 
Sansm™' m ^ ne re gul ar way with the use of relations (183]) . (184")) and the transformation 
matrix connecting the set {F^ sm } (I77|) with the corresponding new set. 



4-3. The asymptotics of radial functions in terms of S, 



a'n's'm' 
ansm 



According to (jHED, ([BHD the wave function ^c^m defined by the asymptotical condition 
(|82|) can be presented in the form of the HSCS-decomposition 

Vansm = ^(l + £ fSLip^ip^l (85) 

j'm' 

where the radial functions are are given by 

fin'smiP) = J ^W(P|^)P 5/2 (1 + t 2 )- 3 ^ ansm (p\Q)dQ. (86) 

Using this relation and the asymptotical formula (|82p for ^ ansm we can express the 
asymptotics of /^ m (p) in terms of the matrix elements S"' 1 ^^ m ' . To this end we have 
to rewrite the expression (I8"2"j) in the hyperspherical coordinates using the relation 

R a = (2M a )- 1 / 2 p + 0(p- 1 ) (87) 

which is valid in the asymptotical region f!72|) . As a result we have 



^ansm(p\ty 



F 2 1/4 M 3/4 



exp [-i(q an p - 7« n )] 



mism i i "a i /o 

Han p 



ST qa'n's'm'j? ol/4 , f 3/4 exp [i{q a > n >p - T^j] 

/ , °ansm r a'n's'm'^ Ivl a ' 1 / 2 



where 



q an = y E — E an , (89) 

7cm = 7cm (/?) = (^a - l)^3- a M a ^ log 2<7 an p. (90) 
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The substitution of anS m{p\^) by the asymptotic expression (1551) in equation 
and the integration with taking into consideration (1761) and (1751) lead to the asymptotical 
formulae for radial functions fansm" 1 ' (p) "where {ansm} is the index of the input channel. 
For indices {jm'} = {a'n's'm'} corresponding to open channels we obtain 

ftnsi m ' (P) = 2(M/i) 3/4 g- n 1/2 {S aa >6 nn >6 aa >6 mm > exp[-i(q an p - % n )} 

p— ►OO 

- (q*n/q«>n>) 1/2 S:' n n ^ m ' exp[i(g QV p - j m )]} . (91) 

For {j, m'} corresponding to closed channels and the continuous spectrum of the 
Hamiltonian h (137|) the result is 

/£ m (p)|^ = 0. (92) 

4-4- The coordination with the boundary conditions of the scattering problem 

The relations (l9"Ij) and ( 1921) show that the radial function fanlm™' {p) does n °t vanish at 
p — > oo only if its indices a'n's'm' correspond to one of open channels. In that case its 
asymptotics contains the single matrix element S^^ m ' . In other words, the function 
describing the motion of the non-interacting fragments in the definite open channel 
(the definite term of the expression ( 1521) ) is the asymptotics of the single term of the 
decomposition (155]) . Other terms of (1551) do not make contribution to that function. 
Thus, the HSCS expansion (1551) is coordinated with the boundary conditions of the 
scattering problem. 

That is not right for the traditional BO approach where R is an adiabatic parameter. 
The reason is that the BO basis functions form an orthogonal system on the surface 
R = const, so that surface appears in the integral for the radial function which replaces 
( l86l) in the BO case. For the calculation of that integral one has to to use for R Q the 
expansion in powers of R~ l . But that expansion, in contrast to the expansion in powers 
of p" 1 (1571) . contains the non- vanishing term depending on the atomic coordinate r a : 

R a = R+ (-lfmsK + m^ 1 ^- R + O (iT 1 ) . (93) 

Consequently, the asymptotic formula for the ingoing (outgoing) wave exp(±ik an R a ) 
contains the factor 



exp 



±ik an m 3 (m a + m 3 ) R r a R . (94) 



The matrix elements of that factor differs from zero for all states of the atom a including 
the excited states corresponding to closed channels. So the radial function does not 
vanish at R —>■ oo for all states with a' = a. It means that all terms of the traditional 
BO decomposition with given a make contribution to the wave functions of all open 
channels with that a. 
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In the numerical calculations it is suitable to use instead of the complex functions $ msm 
( |85l) . ( 1821) their real linear combinations G ansm which contain sm(k an R a — ^an — vrJ/2) 
in the input channel and cos(k a r n /R a r — 7 a / n > — ttJ/2) in all open channels. The HSCS 
representation for G ansm has the form 



Gansm = P^l + ^ £ 9£Ln(p)*i™? . 



(95) 



where the asymptotics of the radial function g a ^ sm (p) with indices {jm'} = {a'n's'm'} 
corresponding to open channel is given by 



'(p) 



p— >oo 



$ aa '$nn'8 ss '8 m m' sm(q an p - ^ an - nJ/2) 



(q an /qo 



Ki 



cos(q a > n rp - 7 



7T 



J/2). 



For {j, m'} corresponding to closed channels and the continuous spectrum we have 



yansm 



(p) 



0. 



p—>oc 



The relation between if^ m ® and S^™ (JH]) can be written as follows 



■IYS 



l + (-l) J S 



fi 



(96) 
(97) 

(98) 
(99) 

(100) 
(101) 



K = 

S = (-l) J [l-iK]- 1 [l + iK], 

where S and K are the matrices with matrix elements 

Sfi = S{ , Kfi = K{ , f = {a'n's'm 1 }, i = {ansm}. 

The relation between functions ^ an sm and G ansm is given by 

Gi = 2-\Mp)-^q] 12 exp[z7r(J + l)/2] £ [l + (-1) J S 

/ 

5. The scattering problem for the radial system 

For the calculation of if-matrix it necessary to solve the system of equations for the 
radial functions which follows from the Schrodinger equation. To obtain that system in 
the simplest form it is suitable to present the AHS Hamiltonian ft, ad (135]) as the sum 

h ad = p~ 2 (l + t 2 )h(l + t 2 ) + W, (102) 

where the additional potential W(p, t) has no singularities in the coalescence points. It 
is given by the formulae 



W = V + 



'l + t 



2\2 



t-ti 



+ 



w (12) +w (13) +w (23) 



(103) 
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w 



.12 



Z 1 Z 2 {2M) 1/2 (1 + t 



,2n1/2 



w 



w 



(13) = z^fi^i + t^ 

\t-tx\ 

(23) = z^y/^i+eyi* 

\t-t 2 \ 



l+i 
1 + 



2\ 3/2 



l+t_ 



2\ 3/2 



3Z 1A /20i-Af 13 ) 



3Z 2 J2(// - M 



2;-w 



(104) 



The three-body Hamiltonian H (1341) takes the form 



d 2 1 

+ - 

P 2 



dp 



;i + t 2 )/i(i + 1 2 ) + (i + t 2 )(j 2 - 2j • i) + 



+ w. 



(105) 



The substitution of ty JKX by the decomposition (1681) in the Schrodinger equation (132]) . 
(I105[) and the projecting onto the state $^ A (p|f2) lead to the system of the differential 
equations for radial functions 



dp 2 4p 2 



j 

i fjm—l( n \ I rp fjm+l( \ 

im,jm—lj \P ) ' -L jm+l,imj \P) 



f d dO m R m \ 



0. 



(106) 



where the sum by j includes both the discrete spectrum and the continuum (see (158]) ). 
The indices J, K, A corresponding to the integrals of motion are omitted. 

The matrix elements P t f, Q%, R%, T imjW , Ug, Wg in (1T06|) are the two- 
dimensional integrals containing the factorized solutions of the CTC problem <f)j m (p\£, fj) 

9 , ,9 



op op 



J vm | Vjm / j 



QW = -QW = (§~ p 

KM = R fM = - 2m 2 ](0 m |(l + t 2 )|0 rm ), 

T imjm ip) = [JfJ+lj-mmf^l + U^x 
x (0 im |(l + t 2 )(-d/d$ + m'cot^l^w), 



(107) 



,J(l+t 2 )/i m (l + t 2 )|0 im ), 



P 

the operator h m is given by the expression 

27T 
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It is evident that the real radial functions gi m ' (p) (1551) also satisfy the system (11061) 
as well as the coefficients of the HSCS decomposition of any solution of the Schrodinger 
equation ( |32i) . 

We did not study the asymptotics of the matrix elements at large p in details as 
it is not necessary for the statement of the scattering problem. However, the principal 
term of the asymptotics of the diagonal matrix element U™ is obtained without any 
complicated calculations: 

U%{p) =E an , i = {ans}. (109) 

p — >oo 

Other matrix elements in (11061) vanish at p —>■ oo. 

To find the matrix K^m™ (ESI) one nas to solve the system (I106P for the real 
radial functions g^ sm {p) for all sets {ansm} (the input channel) corresponding to open 
channels. The solutions have to satisfy the boundary conditions at p = 

flSSm(O) = 0, (110) 



which follow from the finiteness of the CTB wave function G ( 1951) . and the asymptotical 
conditions ( 1961) . (1971) at p — ► oo. The matrix elements S^^ m ' are expressed in terms 
of Karism m ' ^ n accordance with (1991) . The matrix elements of S-matrix can be expressed 
in terms of S^nsm" 1 ' with the help of relation ( I831 . 



6. Concluding remarks 

The suggested HSCS representation is coordinated with the boundary conditions of the 
scattering problem like the AHS one but the HSCS basis functions are essentially simpler. 
The combination of the simplicity and the coordination mentioned is achieved at the 
price of the deviation from the adiabatic idea and the corresponding complication of the 
radial equations. Indeed, the CTC Hamiltonian h (1571) generating our basis functions 
differs from the AHS Hamiltonian h ad (1361) by both the potential energy and the kinetic 
one. Therefore the radial system ( 1106}) contains the non-diagonal matrices U™(p) 
and W™(p) ( 1107!) instead of the diagonal matrix of h ad in the AHS approach. That 
complication, however, looks insignificant as compared with the advantages resulting 
from the properties of the HSCS expansion. First, all calculations are simplified. Second, 
the avoided crossings of the AHS terms are changed by the exact crossings of the CTC 
terms. Third, the use of the well-known and relatively simple CSF allows to perform 
the precise calculations in a wide range of p including extremely large values which are 
essential for the calculation of S'-matrix. These theoretical reasons need, of course, the 
practical verification. 

The essential difference between HSCS and AHS expansions is that the first one 
includes both discrete and continuous spectra while the second one has a pure discrete 
spectrum. In the AHS approach (as well as in the approach of papers [30] - [M]) every 
basis function at large p is localized in the vicinity of the coalescent point where it 



Hyperspherical Coulomb spheroidal representation in the three-body problem 



19 



coincides with the wave function of the bound state of the corresponding atom. Thus, 
that basis is not suitable for the representation of the CTB wave function at energies 
above the threshold of the three-particle breakup. In contrast to the AHS basis, the 
HSCS one includes the functions which are not localized in the vicinities of the coalescent 
points at any p. That allows to hope that the HSCS basis can be used for the adequate 
representation of the CTB wave function at energies above the threshold mentioned. 

The method suggested, by analogy with the adiabatic hyperspherical approach, can 
be named "the diabatic hyperspherical approach" as the avoided crossings of the AHS 
terms correspond to the exact crossings of the CTC terms. In this connection it is 
interesting to develop the semiclassical version of the HSCS approach. 

It is interesting also to apply the HSCS expansion to the calculation of the bound 
states of the CTB systems. That problem is simpler than the scattering one as all radial 
functions vanish at p — > oo. 

The applied three-pole coordinate system on the sphere p = const which is 
connected with the elliptic coordinate system on the plane by means of the stereographic 
projection may be interesting for other three-body problems in the hyperspherical 
approach as it presents the natural way to take into consideration three coalescent 
points of the three-body system. 
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